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ABSTRACT-The results of large scale ab Initio calculations of the rates for 
production of 0 ( 1 D) by dissociative recombination of Oj are presented for 
electron temperatures in the range 100-3000K. A *A U state is the dominant 
dissociative route from v=0 and a 3 Iu state is the most important route from 
v=1 and v=2 . The calculated total rate for 0( 1 D) production from v=0 is 
2.21 ( + 0.21 > -0.24)x10 _7 x(T e /300) - * 1,s near room temperature. The v=1 and v=2 
rates are about 17$ and 47$ smaller respectively, than the v=0 rate at 300K. 
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I. Introduction 


For the molecular ion, AB + , dissociative recombination (DR) with an electron 
is described by 


AB + + e -» A + B (1 ) 

where A or B are ground or excited state atoms or molecular fragments. The 
study of the process of DR of molecular ions with electrons began with the 
suggestion of Bates and Massey (1947) that DR could be an important sink for 
electrons in the ionosphere. It is now recognized that DR can be a rapid 
process in a wide variety of plasma environments. 

Alex Dalgarno has made many significant contributions to the study of DR. 
These have ranged from fundamental studies of the DR mechanism to applications 
in interstellar and atmospheric chemistry. The interstellar studies (Dalgarno, 
1976, 1982, 1983, 1985, 1986; Dalgarno et al., 1973; Black and Dalgarno, 1977; 
Dalgarno and Black, 1976) have explored the DR of a wide variety of molecules, 
ranging in size from H 2 to C 2 H s 0H 2 + . Several of the atmospheric studies of DR 
(Dalgarno, 1958, 1983) have focused on the role of DR in the production of the 
red line of atomic 0 (Dalgarno, 1958; Dalgarno and Walker, 1964; Abreu et al. , 
1986). Additional contributions relevant to the theory of DR have encompassed 
the calculation of the energies (Dalgarno and Drake, 1971) and widths 
(Bransden and Dalgarno, 1952, 1953, 1953a) of autoionizing states and the 
definition of diabatic states (Dalgarno, 1980). An analysis of the mechanism 
of DR (Bates and Dalgarno, 1962), based upon the earlier work of Bates (1950, 
1950a), provided a number of important insights into DR. Perhaps most relevant 
to the work presented below was the prediction that the electron temperature 
dependence of the DR rate coefficient is Tg° • 5 unless the curve crossing of 
the dissociative and molecular ion states leads to a rapid variation of the 
Franck Condon factors with electron energy. 

One of the examples presented by Bates and Dalgarno (1962) was the DR of 0 2 
with an electron, 0 2 + e - -» 0 + 0. It is now well recognized that DR of 0 2 is 
an important source of 0(*D) in the aurora (Dalgarno and Reid, 1969; Bates, 
1982; Rees and Roble, 1986) and airglow (Schaeffer et al. , 1972; Bates, 1982; 
Sharp, 1986). 
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This paper reports the first large scale ab initio calculations of the rate 
coefficient for generation of O^D) atoms from DR of 0 2 in one of the lowest 
three vibrational levels. The next section has a discussion of the techniques 
used to accurately determine the potential curves for the important 
dissociative routes of 0 2 . The method used for calculating the cross sections 
and rates is discussed in Section III. The approach used for determining 
autoionization widths is summarized in section IV and the resulting rates for 
CK'D) production are presented in Section V. 


II. Potential Curves 

In a seminal paper, Bates (1950a) pointed out that DR rates of significant 
magnitude at low electron energies require that a repulsive potential curve, 
describing the motion of the products in (1), crosses the ion curve within the 
turning points of the ion vibrational level undergoing recombination. Since 
the time for dissociation can often be rapid compared to the time for emission 
of the captured electron, the rate of DR can indeed be fast. The process is 
illustrated in Fig. 1 where an ion in the v=1 vibrational level undergoes DR 
by capture of an electron of energy e. For the DR of 0 2 , inspection of the 
relevant potential curves (Guberman, 1983) shows that for the lowest three 
vibrational levels of the ion, only five repulsive states can yield CK'D). The 
states of interest are the lowest states of 3 Iu> 1 ^u* l H g » an 8 3 ^g 

symmetries. The latter two states do not dissociate adiabatically to *D atoms 
but have avoided crossings with higher states that lead to *D atoms. However 
the two n states have exceedingly small calculated widths, .0030eV for 1 Ilg and 
.0022eV for 3 JI g (Guberman, 1987a) at R=2.28l9ao. These widths are about two 
orders of magnitude smaller than the smallest width calculated (see Section 
IV.) for the other three states. Since these small widths will lead to very 
small DR rates, these states are neglected in the calculation of the total DR 
rate for production of *D atoms. The calculated potential curves for the 
remaining three states are shown in Fig. 1 . The potential curves have been 
calculated using contracted Gaussian basis sets of [6s,3p,2d, 1 f] size. The 
basis sets have been chosen to exclude Rydberg character and the resulting 
potential curves describe the diabatic valence states needed for DR. The 
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orbitals are determined separately for each state in Complete Active Space 
Self-Consistent Field (CASSCF) calculations (Roos et al., 1980) and the 
Configuration Interaction (Cl) wave functions are generated by taking all 
single and double excitations from a multireference set of configurations to 
the full virtual orbital space. The direct Cl method (Siegbahn, 1980) and the 
MOLECULE (Almlof)- SWEDEN (Siegbahn et al.) programs were used for the 
calculation of the potential curves. The Cl energies have been corrected for 
the missing quadruple excitations (Langhoff and Davidson, 1974) . The final Cl 
wave functions for the 3 Iy, *A U , and *^u states had 228,036, 185, A00 and 
1 39 ,946 terms respectively. Further details of the calculations will be 
reported separately (Guberman, 1987b). 

The 3 Zy state shown in Fig. 1 intersects the ion at the small R turning 
point of the v=0 level and is the diabatie continuation of the repulsive wall 
of the upper state of the well known Schumann-Runge bands. Dissociation along 
3 Zu after capture by a zero energy electron from v=0 leads to a ground state 
atom and a *D atom, each with 2.5 eV kinetic energy. Dissociation along the 
‘Ay route from v=0 leads to two *D atoms, each with 1.5eV kinetic energy. 
Dissociation along 1 leads to a l S and a l D atom, each with 0.A eV kinetic 
energy from v=0. 


III. Methods for Cross Sections and Rates 

The cross sections, a v , for direct DR (Bardsley, 1968) from ion vibrational 
level v have been calculated using the expression derived by Giusti (1980) 
from quantum defect theory for the case of a single dissociative state and 
many open ion vibrational levels, v' , 

o v =(TTr/2k§)[45 v /(1+£ v ,S v ') 2 ] (2) 


where k| = 2me, k e is the wave number of the incident electron, m is the 
electron mass, and r the ratio of multiplicities of the neutral and ion 
states. The matrix elements £ v are given by 



S v =(TT/2)[<X v |r 1/2 (R) | x d >] 2 


(3) 


where X v and X^ are vibrational wave functions for the ion and dissociative 
states, respectively. The integration in (3) is over both R and the electronic 
coordinates. The vibrational wave functions are determined by the Numerov 
method on a O.OOIao grid between 0.0 and 8.0ao. The integration over R in (3) 
is done by Simpson's rule. The width, T(R), connects the electronic continuum 
of the free electron and ion to the repulsive dissociating resonance state and 
is given by the Golden Rule, 

r(R) = 2TT<{$ ion (x,R)<t> £ (x,R)} |H|$ d (x,R)> 2 . (4) 

In (4) H is the Hamiltonian, the integration is over the electronic 
coordinates represented by x, and the wave functions from left to right in (4) 
are for the ion core, the free electron, and the dissociative neutral state. 
The antisymmetrized product of the ion core and free electron wave functions 
is denoted by { } . 

The indirect DR process (Bardsley, 1968a; Giusti, 1980; Giusti-Suzor et al., 
1983) in which recombination occurs into a vibrationally excited Rydberg state 
which is predissociated by the repulsive routes discussed here, has not been 
included in the cross section calculations. The indirect process is not 
expected to be a significant contributor to O^D) production. 

Rate constants have been calculated by averaging the calculated cross 
sections over a Maxwellian distribution of electron energies. 

IV. Electronic Widths 

The widths given in Eq. (4) have been calculated by using high Rydberg 
orbitals to represent the free electron (Guberman, 1986, 1987a). Since all the 
orbital amplitudes in (4), except for that of the free electron, fall off 
rapidly with increasing x, only the free electron amplitude near the nuclei is 
needed for the width calculation. Except for a normalization factor, high 
Rydberg orbitals have a shape similar to free electron orbitals near the 
nuclei. A density of states is then inserted into (4) for the normalization 
(Guberman, 1987a). The Rydberg orbitals are calculated using the Improved 
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Virtual Orbital (IVO) method (Hunt et al. , 1969) in a basis set of 18 diffuse 
Gaussians positioned midway between the nuclei (Guberman, 1987a). Each Rydberg 
orbital is an eigenfunction of an IVO Hamiltonian which is unique for each 
molecular symmetry. 

The wave function of the ion plus a free electron in (4) is represented by 
antisymmetrized products of a valence Cl description of the ground state of Oj 
with the Rydberg orbital and with the valence virtual orbitals of the same 
symmetry as the Rydberg orbital. For the width calculations, the dissociative 
states were represented by full valence Cl wave functions in which excitations 
out of the 2o orbitals were included. For 3 Ey this procedure leads to a 48 
term Cl wave function for the dissociative state and a 954 term Cl wave 
function for the entrance channel. For ‘A u and the entrance and exit 

channel wave functions have 42 and 630 terms respectively. The size of the 
entrance channel wave function is increased further by allowing the entrance 
channel to mix in higher nonphysical routes from the exit channel. A Feshbach 
(1958) projection operator procedure has been used to implement this procedure 
(Guberman, 1986, 1987a). 

The entire procedure is repeated for several values of the Rydberg orbital 
principal quantum number with the highest accurately determined orbital having 
n=8. The widths calculated from each value of n are then extrapolated to the 
continuum. A comparison to experimentally derived (Galluser et al., 1982) 
predissociation matrix elements for NO indicates that the widths calculated by 
this procedure deviate from experiment by about 15?. 

For 3 Zy, the calculated width is 0.81 eV near R=2.08ao and decreases to 
0.54eV at 1.78ao. For J A U and *1^, the width is relatively constant with R, 
having a value of 0.44eV and 0.29eV respectively at R=2.28ao .varying by only 
about 2 % over the region of significant nuclear overlap with the v=0,1,2 ion 
wave functions. Since the estimated uncertainty in the calculated widths is 
15?, the variation of the width with R has been ignored for and 1 but 
has been accounted for in expression (3) for 3 I^ (S. L. Guberman, 1987b). 


V. Calculated DR Rates 

The calculated rate constants for the three dissociative routes discussed 
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here are shown in Fig.s 2, 3, and 4 for each of the lowest three vibrational 
levels of the nonrotating ion. For 200£T e £400K, the DR rates are given in 
Table 1 . The uncertainties in the rates are average uncertainties computed 
from an estimated ±15? uncertainty in the calculated widths and an estimated 
±0.1 eV uncertainty in the position of the potential curves relative to the ion 
curve. The highest DR rates are from 3 Z^. However, for v=0, the rate along *A U 
is more than half the DR rate along 3 Z^. Since *A U leads to two *D atoms, the 
highest rates for producing *D atoms from v=0 are from *A U . However, from v=1 
and v=2 (Fig.s 2-4), the highest rates for 0(*D) production are from 3 Zu* 

As pointed out by Bates and Dalgarno( 1962) , the temperature dependence of 
the direct DR process will vary as Tg° * 5 if the amplitude of the ion wave 
function does not vary substantially over the range of electron energies of 
interest. Such a situation can arise for any ion wave function if the 
dissociative curve is very steep. In general, an intersection of the 
dissociative route with the ion vibrational level near the peak of the ion 
wave function will lead to slowly varying nuclear overlap terms and electron 
temperature dependencies near -0.5, while intersection at a point where the 
absolute value of the ion wave function amplitude is increasing (decreasing) 
will usually lead to temperature dependencies more positive (negative) than - 
0.5. Note that the denominator in (2) has a tempering effect upon these 
conclusions. Two such cases are illustrated in Cases a and b of Fig. 2 of 
Bates and Dalgarno (1962). The latter Case shows a dissociative route 
intersecting the v=0 level in the nonclassical region on the large R side of 
the ion equilibrium internuclear distance and resembles the situation found 
here for ‘zj,. 

For 'Zy, the amplitudes of the ion wave functions at low electron energies 
increase most rapidly for v=0. The increase for v=1 is less rapid than for 
v=0. For v=2, the intersection is near the outer peak in the v=2 wave 
function. The calculated temperature dependencies given in Table 1 for 1 Zy 
show the predicted pattern. The most positive exponent is for v=0 followed by 
v=1 and v=2. For the other two dissociative routes the exponents are near 0.5 
for v=0 since the intersections with the ion are within the classical turning 
points and the amplitude is changing slowly. For v=1 and 2 and 3 Zy and l A u , 
one must take into account the overlap of the nodes in the ion wave functions 
with the nodes in the continuum wave functions. The exponents are difficult to 
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predict in these cases by visual inspection of Fig.1. 

Also shown in Fig.s 2-4 are the rates for the equilibrium case where the 
vibrational temperature is set equal to the electron temperature. For 3 Zu> the 
equilibrium case rates are nearly indistinguishable from the v=0 rates over 
the span of electron temperatures plotted. For the equilibrium rates fall 

below the v=0 rates at high electron temperatures. For *££, the equilibrium 
rates increase dramatically above 450K, since the rates for v=1 and v=2 are 
substantially larger than for v=0 (Guberman, 1987). 

It is interesting to examine the *D production in each channel as a function 
of ion vibrational level. From v=0, 56? of the *D atoms arise from the *A U 
channel and 44? arise from the 3 Z[j channel. For v=1 , 58? of the *D atoms are 
in the 3 Zu channel and 4? appear in the 1 Ey channel with 0.52eV kinetic 
energy. For v=2 f 45? of the atoms are in the 2.7eV 3 Zy channel and 21? are in 
the 0.64eV *Ey channel. Each vibrational level will lead to a unique three 
component Doppler line shape. From these data, the vibrational distribution of 
the recombining ion could be determined in an experiment that accurately 
determines the Doppler widths of the radiating *D atoms. 

The total rate for producing *D atoms for T e in the above range from v=0 
(Table 1) is 2.21 (+0.21 ,-0.24) x 10 -7 x (T e /300) _o • 1 * 6 cm 3 /sec and is in good 
agreement with the value originally reported by Zipf (1970) of 1.9 x 10~ 7 
cm 3 /sec. Note that the total calculated rate decreases with increasing v. The 
rate from v=2 is about half the rate from v=0. The variation of the rate with 
v contrasts sharply to that for producing 0( 1 S) (Guberman, 1987) shown in the 
fourth column of Table 1 and Fig. 4. The rotationless 0( 1 S) rate for v=2 is a 
factor of 80 larger than the rate from v=0. The total calculated rate up to 
T e =3000K for *D production is shown in Fig. 5. Rates for both J=0 and 
rotationally averaged rates in which the rotational temperature is the same as 
T e are shown. Clearly, accounting for the rotational temperature has only a 
slight effect upon the calculated rates. The rotationally averaged and the 
rotationless results from v=0 and for the case of equal vibrational and 
electron temperatures are each superimposed on the scale of Fig. 5. 

There have been several aeronomic determinations of the *D production rate 
by DR. The literature prior to 1982 has been reviewed by M. R. Torr and D. G. 
Torr (1982). Among the deduced production rates are (2 . 1±0 . 48)x1 0~ 7 x 
(T e /300 )~- 55 (P. B. Hays et al, 1978; R. Link et al. , 1981) and 



(1 .96±0.67)x10 -7 x(T e /300)“* 55 (Abreu et al. , 1983). These rates were obtained 
by combining the total rate found by Walls and Dunn (1974) and Torr et al. 
(1976) with a model determination of the *D quantum yield. A recent 
determination (Abreu et al. , 1986) which improves the models by including 
quenching of O^D) by 0( 3 P) leads to a rate of 1.9x10 7 x(T e /300)~* 5S . The 
deduced rates are in excellent agreement with the v=0 rate calculated here at 
room temperature. The atmospheric quantum yields are derived at ionospheric 
temperatures. At T e =760K, the *D production rate calculated here is 
1 . 41 (+0. 1 4,-0 . 1 6)x1 0~ 7 compared to the model rate (Abreu et al., 1986) of 1.14 
x 10“\ Considering the likely precision on the latter rate, the rates 
reported here are in agreement with the model rates. Note that the precision 
on the rates calculated here do not exclude the possibility that the models 
apply to vibrationally excited 0^. 

None of the above models have included *D production due to cascade from *S. 
The DR contribution to cascade will be unimportant if the ion is mostly in v=0 
but should be considered if a substantial fraction of the ions are 
vibrationally excited. 


VI. Conclusions 

Three neutral dissociative states of 3 Iu, 1 A U , and 1 symmetries have been 
identified as the major sources of 0( 1 D) from DR of the lowest three 
vibrational levels of Oj. DR rates along each of these routes have been 
calculated for electron temperatures in the range 100-3000K. The total rate 
for generating DR from v=0, 2 . 21 (+0. 21 ,-0 . 24) x 10" 7 x (T e /300)~° • 4 6 em 3 /sec, 
agrees with prior rates determined from aeronomic models and laboratory 
measurements. However, the variation of the rate with v and the precision on 
the calculated rates leaves room for the possibility that the atmospheric 
models and the laboratory experiments could be observing vibrationally excited 
ions. The total rate for generating 0(*D) decreases with increasing v in 
contrast to the rates for generating 0( 1 S) which increase markedly with 
increasing v. 
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TABLE 1 . Dissociative recombination rates along three channels for each of the lowest three vibrational levels of the 
ion at 200K<T e £H00K and J=0.* 
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Figure Captions 


1. Potential energy curves for DR leading to O^D). The three important 
calculated repulsive routes for DR of the lowest three vibrational levels are 
shown. The bound potential for the ion is taken from Krupenie (1972) and 
calculated points (Guberman, to be published). The vibrational wave functions 
are each plotted on the same scale and are obtained from the numerical 
solution of the nuclear Schroedinger equation. 

2. Calculated DR rates along the 3 Zu dissociative route leading to 0( 3 P) + 

0( X D) from the lowest three vibrational levels of the ion. Also shown is the 
rate for the case in which the vibrational temperature is the same as the 
electron temperature. The rates are for non-rotating molecules. 

3. Calculated DR rates along the X A U dissociative route leading to 0 ( 1 D ) + 
0(*D) from the lowest three vibrational levels of the ion. Also shown is the 
rate for the case in which the vibrational temperature is the same as the 
electron temperature. The rates are for non-rotating molecules. 

4. Calculated DR rates along the l £y dissociative route leading to 0( X S) + 
0( X D) from the lowest three vibrational levels of the ion. Also shown is the 
rate for the case in which the vibrational temperature is the same as the 
electron temperature. The rates are for non-rotating molecules. 

5. Calculated total DR rates for production of 0(*D) from the lowest three 
vibrational levels of the ion. Both the J=0 rates and rates for which the 
rotational temperature is the same as the electron temperature are shown. 
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